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Among Markovian processes, the hallmark of Levy flights is superdiffusion, or faster-than- 
Brownian dynamics. Here we show that Levy laws, as well as Gaussians, can also be the limit 
distributions of processes with long range memory that exhibit very slow diffusion, logarithmic 
in time. These processes are path-dependent and anomalous motion emerges from frequent re¬ 
locations to already visited sites. We show how the Central Limit Theorem is modified in this 
context, keeping the usual distinction between analytic and non-analytic characteristic functions. A 
fluctuation-dissipation relation is also derived. Our results may have important applications in the 
study of animal and human displacements. 


PACS numbers: 05.40.Fb, 89.75.Fb, 87.23.Ge 

I. INTRODUCTION 

Levy flights (LFs) represent one of the most impor¬ 
tant extensions of the Central Limit Theorem (CLT), a 
cornerstone of probability theory LFs are sums 

of independent and identically distributed random vari¬ 
ables that admit non-Gaussian limit laws due to their 
very large fluctuations. They find physical applications 
in laser cooling Q, optics Q or chaotic transport Q. 
LFs are also paradigmatic of superdiffusive processes, 
i.e., anomalous types of transport where the character¬ 
istic diffusive length scale l{t) of an individual particle 
grows with time as with a > 1/2, that is, faster than 
in the classical Brownian motion (BM) il. 

In recent years, LFs (as well as the related Levy walks 
[l^ ) have become prominent for modeling diffusion in 
a variety of complex systems. Power-law distributions 
of step lengths with diverging variance, a key feature of 
Levy processes, are found to describe well the trajecto¬ 
ries of immune cells in the brain , the displacements 
of animals [l^ - fTsjl and hunter-gatherers [H, [T3] in their 
environments, or the travels of modern humans within 
and between cities [Tsl - f^ . However, the assumption 
of independence between steps does limit the applica¬ 
bility of genuine Levy processes for modeling real sys¬ 
tems, where non-Markovian effects and correlations can 
be strong. Deeper analysis of empirical data actually 
reveals that the diffusion of humans and animals (even 
those exhibiting Levy patterns) is in general sttfediffusive 
at large times, i.e., with [2ll - [25l| . Further¬ 

more, l{t) commonly grows more slowly than a power-law 
of time, namely, in a logarithmic way [^[^,1^: this be¬ 
havior is even in sharper contrast with the superdiffusion 
of simple LFs. 

Logarithmic diffusion can be generated in several ways. 
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for instance, by continuous time random walks mod¬ 
els with superheavy-tailed distributions of waiting times 
[ 2 ^ , or by certain iterated maps [13, . In the context 

of animal and human mobility, an important but little 
explored mechanism that may lead to very slow subdif¬ 
fusion is spatial memory: many living organisms actually 
keep revisiting familiar places [29l . l30j| . Here, we 

seek to understand, with the help of a solvable model, 
how this type of memory can act as a self-attracting force 
which drastically constrains diffusion towards limited ar¬ 
eas, giving rise to “home ranges”, and how this property 
can still be compatible with power-law distributed step 
lengths. 

The dynamics and limit distributions of constrained 
LFs are not well understood, except for processes sub¬ 
jected to long waiting times or in external potentials, 
mainly [3, [3l| . Several limit theorems also exist for spe¬ 
cific problems of sums of correlated random variables 
[^ . and a few random walks with infinite memory of 
their previous displacements have exactly solvable first 
moments [33 - [35l| . Yet, very little is known on LFs com¬ 
posed of non-independent steps, in particular processes 
with self-attraction. Self-attracting random walks are 
path-dependent processes where a walker tends to re¬ 
turn to previously visited sites [H, [33|. Numerical sim¬ 
ulations and scaling arguments clearly show that self- 
attracting walks can exhibit subdiffusion (38l - l40l| . These 
mathematically challenging processes cannot be readily 
analyzed with better known frameworks for subdiffusive 
phenomena, such as fractional Fokker-Planck equations 
Q or scaled Brownian motions [13, . They are more 

related with diffusion in quenched disordered media 
where some rigorous connections have been made with 
the Sinai model [i^ ]. 

In this study, we heuristically modify the CLT for 
processes that exhibit very slow diffusion, and show 
that such modification exactly describes a class of self- 
attracting LF and self-attracting random walks. The 
characteristic functions having a similar structure than 
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in the ordinary CLT, Gaussian and Levy distributions 
emerge asymptotically in space, although the dynamics 
is strongly subdiffusive. We also derive a fluctuation- 
dissipation relation in the Gaussian case. 

II. GENERAL FORMULATION 

Let P{n, t) be the probability that the position Xt of a 
particle at time t is n (where n and t are discrete), given 
that the particle is located at the origin n = 0 at t = 0. 
We consider discrete, one dimensional walks, keeping in 
mind that discreteness is not relevant in the asymptotic 
limit. The results can also be extended higher dimensions 
straightforwardly. 

We recall that for a standard random walk com¬ 
posed of t i.i.d. displacements t'i with distribution p(f), 
the characteristic function of Xt, defined as P{k,t) = 
S^-oo takes the form 

=p(fc)‘(1) 

where p{k) is the characteristic function of £. Since 
p{0) = 1 by normalization, in the unbiased {{£) = 0) 
and symmetric case, an expansion near k = 0 gives: 

p(fc) = 1-C|fcr+ ... (2) 

Two basic situations emerge: the analytic case fi = 2, 
corresponding to (£^) < oo (and C = (£^)/2), and the 
non-analytic case 0 < p < 2 when does not exist, 
due to a power-law decay of p{£): 

p{i) ^ (3) 

at large i Gombining 0-0 yields the celebrated 

Gaussian-Levy CLT: 

P(fc,t)(4) 

Eq. 0 implies a scaling law P{n,t) —>■ t 
where the scaling function f{x) is a Gaussian or a sym¬ 
metric Levy law L^ q(x), for /i = 2 and 0 < /r < 2, re¬ 
spectively. The latter case is superdiffusive as the typical 
diffusion length is oc 

Consider now a simple modification of Eq. ©: sup¬ 
pose that for certain diffusion processes with memory or 
sums of correlated random variables (we do not need to 
specify a model at this point), P is not an exponential 
function of t but a power-law: 

P(fc,t) ~ (5) 

at large t and small k. The function a{k) satisfies 
a(0) = 0, owing to the normalization P{k = 0,t) = 1. 
Again, a{k) can be generically analytic or non-analytic 
near k = 0. In the first case, since P{k,t)* = P{—k,t) 
and \P{k,t)\ < 1, the Taylor expansion of the exponent 
must be of the form a{k) ~ iaik -\- a 2 k'^ ..., with oi and 


02 two real constants and 02 > 0. For simplicity, we first 
consider oi = 0, or motion without bias. 

In the non-analytic case, the same arguments lead to 
a{k) ~ afi\k\^ with 0 < p < 2 a priori, and o^ > 0. 
Inserting into 0, we see that the main difference with 
0 is that the variable t is substituted by Int. Hence: 

where the limit laws f^iix) are the same as in the or¬ 
dinary CLT. If /i = 2, diffusion is Gaussian but very 
slow: {X^) = 202 In t, in sharp contrast with BM, where 
{Xt) = 2Dt. [In this case, Eq. 0 should not be con¬ 
fused with the log-normal distribution, where the log¬ 
arithm applies to the space variable, not the temporal 
one.] A basic Markovian example is, by construction, a 
scaled Brownian motion, which is a BM where the time 
T is rescaled as t = e^. Such process is also equivalent to 
a BM with a time-dependent diffusion coefficient, D{t), 
decaying as Ijt at large t (i^ . 

In the non-analytic case, the situation looks paradoxi¬ 
cal at first sight. The ensemble average (Xf) = 00 like in 
ordinary Levy processes due to the broad tails of 
(or due to the fact that d'^P{k,t)/dk'^ does not exist at 
k = 0, from Eq.0). Yet, Eq. 0 also defines a typical 
diffusion length l{t) oc (Int)^/'^, which grows extremely 
slowly. Therefore, based on this scaling length l{t), mo¬ 
tion is strongly subdiffusive and all the finite moments, 
{\Xt\'') with V < p, also evolve very slowly, as 
Still, the process keeps superdiffusive features through 
the divergence of the second moment. This situation is 
reminiscent of scaling violation, which also arises in con¬ 
tinuous time random walks or Levy walks [l^ . 1^ . 


III. RANDOM WALKS WITH RELOCATIONS 

We now consider a concrete class of non-Markovian 
walks for which the above ideas apply. The processes of 
interest are self-attracting, namely, they tend to revisit 
locations visited in the past. Particular examples were 
studied numerically in 1221 123^ as animal movement mod¬ 
els, or theoretically in [^, 1^. We present here a unified 
view of this class of processes. 

Let g be a parameter (0 < g < 1). At any time t, the 
walker chooses its next position according to the follow¬ 
ing rules: 

(i) with probability 1 — q, it performs a random dis¬ 
placement i drawn from a given distribution p{£) like in 
standard random walks or Levy flights; 

(ii) with the complementary probability q, it jumps (or 
’reset’) directly to the site occupied at some previous time 
t' <t. The time t' is chosen according to a given prob¬ 
ability 7rt(<'), or memory function, with X]t'=o = 1 
by normalization. 

The rules are depicted in Fig. 0, with two simulated 
examples in Fig. 0. Note that in (ii), the next target 
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a) 


b) 




FIG. 1: (Color online) a) Schematic view of a process re¬ 
locating at a constant rate (q) to sites occupied at previous 
times, these times being chosen stochastically. The numbers 
label the beginning and end of each excursion. Each end is 
followed by the beginning of the next excursion (arrow), b) 
Two simulated trajectories corresponding to Levy excursions 
with p{£) , relocation rate q — 0.05 and relocation 

kernel given by Eq. ([HI) [panels at the same scale]. 


site is chosen independently of its distance to the location 
Xt of the walker. If Trt{t') = St'fi, the site chosen for 
revisit is unique (the origin), a case which corresponds to 
the well-studied random walk with resetting to the origin 
[ril - IHflj l. For more general kernels, the walk is strongly 
path-dependent but still described by a master equation: 

OO t 

P{n,t+1) = (l-q) p{i)P{n-e,t)+qY'^t{t')P{n,t'). 

£— — oo t'—O 

(7) 

Standard random walks or Levy flights are recovered for 
q = 0. If g 7 ^ 0, the last term indicates that site n can 
be chosen to be occupied at time t -I- 1, provided it was 
visited at the earlier time t'. 

We first consider a uniform memory function, that is, 
independent of t'\ 

^.(0 = ( 8 ) 

We call this case the preferential visit model (PVM): with 
such kernel, rule (ii) is simply equivalent to choosing a 
given site n (among all visited sites) with probability 
proportional to the number of visits received by n since 
t = 0. Therefore the walker is prone to revisit familiar 
sites, at the expanse of rarely visited ones. The moments 
(Xf^) where calculated in [23 for the PVM with nearest 
neighbor (n.n.) steps {£i = ±1) in rule (i). To solve Eq. 

© more generally, we define the Laplace transform of 
P{k,t): 

OO OO 

P(fc,A) = ^A* ^ e-*'="P(n,t). (9) 

t—0 n——oo 


By taking the double transform of Eq. 0 with the kernel 
dSl) and writing X*/(t -b 1) = A“^ vfdu, we obtain: 

P{k,\) — l = {l—q)p{k)\P{k,\)+q f du ^ (10) 

Jo 1 - w 

Taking the derivative of Eq. (nni), one obtains a first- 
order ODE in the variable A. As P{n,t = 0) = i5o,„, the 
condition P{k,0) = 1 must be enforced, leading to the 
exact solution: 

P{k,X) = (1 - A)-b-“W] [1 - (1 - g)p(fc)A]““^'=^ (11) 
with 


a{k) = (1 - g) 


1 -.P(fc) 

1 - (1 - g)p(fc)' 


( 12 ) 


We can infer the large t behavior of P{k,t) by studying 
the divergence of P{k,X) near A = 1, with k fixed but 
small. Noting that a{k) <C 1, Eq. (fTTl) yields P{k,X) ~ 
(1 — . This expression is simply inverted as: 


P(fc,t) ~ (13) 


as announced in ([5]). In the absence of bias, one can use 
Eq. ([5]), which, combined with (IT^ . gives the exponent: 


a{k) ~ (14) 

q 

implying the limit law ([6]). We conclude that this ran¬ 
dom walk always diffuses logarithmically, unlike other re¬ 
inforced walks that exhibit transitions to localized states 
[ 3 ^ 1^ . Numerical simulations confirm the very slow 
dynamics, even for /r < 2: a perfect agreement with the 
prediction {\Xt\'') ^ (Int)'^/'^ ior v < p, is observed in 
Eig. Importantly, the scaling function f{x) in this 
non-Markovian process is the same as for the underlying 
Markovian process between relocations (or with g = 0). 
This property stems from the fact that the cumulant 
characteristic function Inp(fc) [Eq. (|T])] and the function 
a{k) [Eg. dT^ ] have the same leading behavior at small 
fc, except for a multiplicative constant. In other words, 
the analyticity or non-analyticity of P{k, t) is preserved 
when g is set different from zero. 


IV. GENERALIZATIONS 

We now show that several extensions of the PVM also 
admit a propagator of the form given by Eq. ©• 


A. Decaying memory 

The results of the previous Section do not change quali¬ 
tatively by considering memory kernels other than a pure 
preferential one. Eor instance, the time in the past t' may 
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a) b) 



FIG. 2: (Color online) Preferential visit model in Id. a) 
, obtained from simulations with different ^ and 
u (averages over 5 x 10® runs), is proportional to Inf as ex¬ 
pected. b) Mean and variance of Xt for a n.n. walk with 
bias a in rule (i). Colored solid lines are simulations and 
dark dashed lines, theory, c) Normal diffusion for spatially 
uniform relocations. 


be chosen not uniformly like in Eq. ([5]) but with a proba¬ 
bility decaying with t — t', the interval of time between a 
remembered occupation and the present time. Consider, 
for instance, a power-law memory decay: 


TTtit') 




(15) 


with /3 > 0 an exponent. Here, the visits are still pref¬ 
erential, but with a tendency towards more recent sites 
(an effect actually observed in human mobility [5l|). If 
/3 < 1 the sum in (ITSl) diverges at large t and can be 
substituted by an integral; by taking the Fourier trans¬ 
form of 0 and making the ansatz P{k,t) ~ t one 
obtains an integral equation for a(fc): 


1 - (1 - q)p{k) = q{l - (3) 



y)-/3y-a(fc)^ (16) 


Combining Eqs. HID and 0 gives, at small fc: 


a{k) ^ .F(/3)C|fcr, 


with J-iP) = 


(1-/3) [\u{l 
Jo 


(17) 

-1 

u)“^ln(l/M) 


Eq. (fTTp shows that the scaling law ([ID applies to more 
general processes than the PVM. [Ea. lfTHl is recovered for 
P — 0.] Interestingly, J^(l) = oo, which indicates that 
the scaling form ([SD breaks down for /3 > 1. Actually, 
a similar calculation to the one above shows that, for 
/3 > 2, memory decays too fast to be relevant and the 


usual CLT dlD is recovered. Of course, these results do 
not mean that the aforementioned preservation property 
holds for arbitrary irtit'). For instance, for memory walks 
with 1 < /3 < 2 and steps £i of finite variance, the process 
is non-Gaussian Likewise, Brownian random walks 
and Levy flights subjected to stochastic reseting to the 
origin have asymptotic probability densities which are 
non-Gaussian [47j | and non-Leyy [50j, respectiyely. 


B. Model with bias 


We now study the response of the non-Markoyian 
walks (at fixed q) to the presence of a constant forcing, 
namely, a bias a = {£) = 0. Here, we 

assume < oo or /j, = 2. By taking the first mo¬ 
ment of Eq. ©, an equation for the ayerage position 
(Xt) = obtained: 


(A/t-i-i) — (1 — q)[{Xt) -I- a] + g^ TTt{t'){Xt'), (18) 

t=o 


for any kernel TTt{t'). We now denote as the mean 

square displacement of the walker at zero bias. It is 
easy to show that {Xf)g obeys exactly the same equa¬ 
tion as dUD, where a has to be replaced by {£‘^)g = 
with Pg{£) unbiased. We deduce an Ein¬ 
stein fluctuation-dissipation relation (FDR): 

(W) = J^jXf)g (19) 

The exact equality (ITOl) is general: it is yalid at all t 
and for any kernel T^t{t') (allowing to recoyer results on 
the resetting to the origin with bias [3). Despite of 
being out-of-equilibrium, the FDR with constant bias in 
this system is the same as for ordinary random walks, 
where the response {Xt) is entirely determined by the 
fluctuations at zero bias. With the kernel (ITSl) and /3 < 1, 
the drift is thus logarithmic: {Xt) ~ aX^iF{p) Int, from 
Eqs. m and (ITTP with p = 2. The time evolution of the 
first moment {Xt) is displayed in Figl^Jo-left for different 
parameter values. 

In other words, the effective friction coefficient of the 
walker (oc a{Xt)~^) grows linearly with t. This illus¬ 
trates the non-stationarity emerging from long range 
memory and the increasingly sluggish dynamics caused 
by frequent relocations to the same preferred sites. 

We further show that the combination of memory and 
bias has a drastic impact on the fluctuations of Xt around 
{Xt). We take, for example, the PVM with n.n. steps 
in rule (i), and expand Eq. (1121) . which is valid for any 
p{£), near k = 0. Now using p{k) = 1 — iak — ^k^ -|- ... we 

obtain P{k,t) ~ exp[—f/ijA: — ^atk^], which corresponds 
for P{n, t) to a Gaussian of mean pt and variance at. We 
recover pt = Int, see (IT^ . and obtain for ap. 





Inf. 


( 20 ) 
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If q is small, the presence of a bias therefore strongly 
amplifies the fluctuations of Xt, as the 2"'^ term in (1^01) 
is > 0 and dominant. This effect is displayed in Fig. [2b- 
right. For ordinary n.n. random walks, on the contrary, 
the bias decreases the fluctuations: in that case at = 
(1 — a^)t and motion becomes deterministic at a = 1 
(see e.g. [s^l. 

V. DISCUSSION AND CONCLUSION 

In summary, we have shown that Levy and Gaussian 
distributions can emerge generically far from the domain 
of applicability of the CLT, namely, in strongly subd- 
iffusive path-dependent processes. We emphasize that 
the processes studied here exhibit subdiffusion because 
the relocation sites are selected heterogeneously in space. 
This situation is also encountered in the reseting to the 
origin, an extreme case where only one site receives all re¬ 
locations, causing the typical diffusion length l{t) to tend 
to a constant . To illustrate the importance of uneven 
relocations, one may by contrast consider a n.n. random 
walk, which, in rule (ii) above, relocates to a site cho¬ 
sen randomly and uniformly among the visited sites. In 
this case, l{t) roughly obeys dl/dt ^ (2i?/Z)[(i?/2)/(I/g)], 
with R = \J2D jq the characteristic diffusion scale be¬ 
tween two relocations, 2i?/Z being the probability of re¬ 
seting near the edges of the territory covered by the walk. 
This leads to lit') ^ y/ADt, a normal diffusive behavior, 
which is qualitatively confirmed by the numerical simu¬ 
lations of Figure HJ:. 

The emergence of logarithmic diffusion can be under¬ 
stood qualitatively by drawing, from Fig. [Tb, an anal¬ 
ogy with a branching random walk (see, e.g. [13|)- 
Consider an initial normal random walk with a constant 
branching rate qb. At each branching event, a new ran¬ 
dom walk is created which starts from the current posi¬ 
tion of the parent walk. The walks are independent, do 
not disappear, and all branch at the same rate qt. The 
process follows until it is stopped at some final time T. 
Let then imagine a single walker starting at the origin 
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